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Abstract 

We present arguments that the structure of the spectrum of the supersymmetric 
matrix model with 16 real supercharges in the large N limit is rather nontrivial, 
involving besides the natural energy scale ~ A 1//3 = (g 2 N) 1 ^ 3 also a lower scale 
~ A^iV -5 / 9 . This allows one to understand a nontrivial behaviour of the mean 
internal energy of the system E oc T 14 / 5 predicted by AdS duality arguments. 



1 Introduction 

The AdS/CFT duality is an efficient method allowing one to obtain nontrivial predictions 
for many observables in certain supersymmetric field theories at strong coupling [1]. Most 
results were derived for M = 4 4D supersymmetric SU (N) Yang-Mills theory in the 't 
Hooft limit N — > oo with fixed and large A = g 2 N. The wonderful Maldacena conjecture 
that the properties of this theory at A 3> 1 can be derived by studying classical solutions 
of 10.D supergravity is not proven now. However, it was verified in several nontrivial 
cases where exact solution is known. Arguably, the most lucid example is the circular 
supersymmetric Wilson loop [2]. For large N, its vacuum average can be perturbatively 
evaluated in any order in A. The sum of the perturbative series is 

w«. = ^ • (i-i) 

On the other hand, the same quantity can be calculated at large A on the AdS side. The 
result 



^circle = \lz^~ 



A3/4 



8^ 
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coincides exactly with the asymptotics and preasymptotics of (II. ip . Another important 
example is the so called cusp anomalous dimension [3]. 

Besides vacuum averages and the scattering amplitudes, one can also calculate ther- 
modynamic characteristics. Thus, the mean energy density of M = 4 AD SYM system 
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at nonzero temperature was evaluated at strong coupling at leading [4] and subleading [5] 
order. The result is 



7T 2 7V 2 A3 45C(3) 
1 ' 4 + 32A^ + 



;i.3) 



(the coefficient in front of T 4 is the coefficient in the Stefan-Boltzmann law). In this case 
the exact result for the function /(A) multiplying the factor 7r 2 iV 2 T 4 /2 in the expression 
for E is not known (though asymptotic expansion (11.31) matches perfectly the known 
perturbative expansion of /(A) at small A), and it is not thus clear how the nontrivial 
coefficient 3/4 in the strong coupling asymptotics is obtained, if staying in the framework 
of field theory and not going to the AdS side. 

Duality relationships can be established and duality predictions can be made, however, 
not only for AD theory, but also for its low-dimensional "sisters" obtained by dimensional 
reduction. In particular, by studying a certain charged black hole solution in 10D su- 
pergravity, one can evaluate the average internal energy of the supersymmetric quantum 
mechanical system obtained from Af = 4 AD SYM by keeping there only zero spatial 
field harmonics. The model involves 8 complex or 16 real supercharges. The M = 4 AD 
SYM model can in turn be obtained by dimensional reduction from TV = 1 10-dimensional 
theory. To distinguish it from the models obtained by dimensional reduction from AD 
and 6D M = 1 theories, we will refer to it as "10-D SQM model". Its Hamiltonian is 



H = X -E«Et + 9 ^j^A\A)A\A] + |rX(r,) a/3 A^ , (1-4) 



where i, j — 1, . . . , 9, a = 1, . . . , N 2 — 1 , and a, (3 — 1, . . . , 16. Ef are canonical momenta 
for the bosonic dynamic variables A". Now, A^ are Majorana fermion variables lying in 
the 16- plets of 50(9). Tj are 9-dimensional (real and symmetric) T-matrices. One can 
introduce 8(iV 2 — 1) holomorphic fermion variables, 

^ = A? - iAS , ^ = K-iK &) (1.5) 

such that the wave functions depend on Af and on [i\ 8 , while fl^ 1 '-* 8 are the fermion 
canonical momenta, Ji = d/dfi. Only the states \1/ satisfying the Gauss law constraint 

G a S> = f abc (a*E? - l -X b a Xlj V = . (1.6) 

should be kept in the spectrum. 

In QD and AD theories, holomorphic fermion variables are defined more naturally as 
Weyl fermions lying in the complex representations of the rotational groups SO (3) and 
50(5) (the spinor representation is real in 50(9) ). For example, in AD theory, the third 
term in the Hamiltonian is 

-^r^A b >,) a % a, = 1, 2 (1.7) 
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(A aQ = d/d\°)- We see that one can define in this case (and also in the 6D case) the 
fermion charge F = X^X aa that commutes with the Hamiltonian. In the 10D case, the 
fermion charge ^ a a Ji aOL is not conserved. 

The coupling constant carries dimension here, [g 2 ] = m 3 . The natural energy scale of 
the theory is thus 

Schax ~ (^iV)V3 = A V3 . (L8) 

The duality prediction for the average internal energy is [6, 7] 



T 
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A question arises whether this rather nontrivial critical behaviour can be understood 
in terms of the dynamics of the system (jl.4p without going to the supergravity side. Even 
though the system (jl.4p is complicated, it is just a QM system with large, but finite (for 
finite N) number of degrees of freedom. The analysis of its dynamics at strong coupling 
is a priori a much more simple task than the analysis of a strongly coupled field theory. 
And, indeed, in recent papers [7,8] (see also [9]) a numerical analysis of the system (11.41) 
was performed. The results are in a good agreement with (jl.9p . Can one understand 
it also analytically (staying firmly on the SQM side) ? Our answer to this question is 
positive. 

However, before giving this answer (it will be presented by the end of the next section), 
we are in a position to describe a proper context where the question should be posed and 
remind some well-known facts. 

2 Thermodynamics and the spectrum 

As a warm-up, consider the harmonic oscillator, H = (p 2 + u 2 x 2 )/2, at finite temperature. 
The partition function is 

At large temperatures T = /3 _1 , 



T 

The latter result can also be obtained semiclassically 



(2.2) 



v high T ~ ^semicl 



/^M> +uV) } = §- (2 - 3) 



1 We sketch the derivation of the leading asymptotics oc T 14 / 5 in the Appendix. 
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The mean energy is 



(E): 



At low temperatures, (E)t ~ uj/2 + 0(e while at high temperatures 



2tanh^ 



(E) 



T . 



(2.4) 



(2.5) 



The behaviour (E) high T oc T is characteristic not only for the oscillator, but for 
any reasonable QM system. Basically, it is the analog of the Stefan-Boltzmann law in 
zero spatial dimensions. For the oscillator with several (many) degrees of freedom this 
number multiplies T in the high-temperature estimate for {E)t- 



2.1 YM matrix models. 

For purely bosonic matrix models with the Hamiltonian being the sum of two first terms in 
(jl.4p . the pattern of the spectrum and the temperature dependence of the average energy 
and other thermodynamic functions is clearly understood and is much simpler than that 
for supersymmetric models. It does not depend much on whether the model is obtained 
by reduction from 4D YM theory (i = 1, 2, 3), 6D theory {i = 1, . . . , 5), or 10D theory 

(' 1 

Let us first estimate the energy of the ground state. The simplest variational Ansatz 

is 

oc exp {-a(A*) 2 } = . (2.6) 

The contribution of the kinetic term in the Hamiltonian to E var is estimated as 

E% « aN 2 , (2.7) 

where we have kept only the dependence on iV > 1, not worrying about the dependence 
on D and about numerical coefficients. (Note that (^4 2 )o ~ N 2 /a in this limit.) To 
estimate the contribution of the potential part, use 



(A a A b ) ~ /u 9 , (A a A"A c A d ) 



var Ar 2 V^- u ; 



N 2 ' ' iV 4 

(irrelevant spatial indices are suppressed). We obtain 

9 2 (A\ ^ g^N 3 
N ~ a 2 

Adding this to (12. 7p and minimizing over a, we find 

a ~ A 1/3 , E ~ iV 2 A 1/3 . (2.9) 

In other words, the estimate for the vacuum energy is obtained by multiplying the natural 
energy scale fll.Sp by the number of degrees of freedom ~ N 2 . The characteristic size of 
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the wave function (" extent of space" in the terminology of Refs. [7,8]) is also determined 
by this scale, 



extent of space ~ - - 2° ~ — ~ — — . (2.10) 

P N 2 a E char 1 ; 

The gap between the first excited state and the vacuum state is also of order -E c har- 

The partition function at high temperatures T ^> A 1 / 3 is easily evaluated by semiclas- 
sical methods [10]. We have 



-0H 



II ^^e"^ n ^ eMiA a Q f abc Atm . (2.11) 

ai a 

The integral is saturated by the characteristic values 



and is estimated as 



J semicl 



T v §7V 2 (D-2) 

AV5 • < 212 > 



which gives [11] 



(£)3t ~ \n\D-2)T . (2.13) 



One can also evaluate corrections to this leading order semiclassical result. To esti- 
mate the next-to-leading correction, one should roughly speaking insert the factor ~ 
f3 2 d 2 Vl{dA a i ) 2 ~ f3 2 X(A a i ) 2 in the integrand in f[2TTTj) (V is the potential) [10,12]. This 
gives □ 



^high T 



A 1/2\13 W2 (M 



^1/3 \ ^3/2 



(2.14) 



We see that the correction is of order one at T ~ A 1 / 3 . At temperatures much less 
than the characteristic spectral gap -E c har ~ A 1//3 , (E)t behaves in the same way as the 
oscillator average energy (12.41) coinciding with the vacuum energy up to exponentially 
small corrections. This pattern was confirmed by numerical calculations [13]. 



2 The parameter c was evaluated numerically for D = 10 in Ref. [11]. It is negative. 
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2.2 Supersymmetric models 

One important feature distinguishing all SYM models (10D, 6D, and AD) from the YM 
models is the presence of continuous spectrum associated with flat directions in the poten- 
tial [14]. The classical potential V ~ g 2 Tr{[Aj, A,] 2 } turns to zero when [Ai, Aj] = 0, 
i.e. when all Ai belong to the Cartan subalgebra. In purely bosonic model, this clas- 
sical degeneracy is lifted by quantum corrections. It supersymmetric models, this does 
not happen. As a result, the states can smear along the flat directions, the motion 
becomes infinite, and the continuous branch of the spectrum exists. For low energies, 
E < A 1/3 , the wave functions of these states can be evaluated in the framework of the 
Born-Oppenheimer approximation [14-19]. In the leading order, they can be chosen in 
the form 

^continuous ~ X(^slow) ^A slow Ofast) ■ (2.15) 

x s i ow in this expression stand for A\ and their fermionic partners. The Cartan subalgebra 
index a runs from 1 to r = N—l. Xf ast are the transverse to the valley components of A and 
their fermionic partners. The motion across the valley is described by the Hamiltonian of 
supersymmetric oscillator with the frequency ~ </|A|. ^A^^fast) is its ground state. The 
function x^siow) is the eigenfunction of the effective Born-Oppenheimer Hamiltonian. To 
leading order, the latter is just the Laplacian 

" B0 ~ + - (2 ' 16) 

such that x ^ exp{ikfA^}. 

In AD and 6D theories where the conserved fermion charge exists, one can ask what is 
its value for the continuum spectrum states (12.151) . Consider for simplicity AD theory. The 
fast oscillator Hamiltonian depends on (D—2)N(N—1) = 2N(N— 1) real bosonic variables 
and 2N(N — 1) holomorphic fermion variables [ iV(iV — 1) being the total number of roots 
in SU(N)]. When N = 2, the ground state of the fast Hamiltonian has the structure [16] 

^(^t) ~ exp |-^l(^ t ) 2 } {X ba X b a + * e bc A b > 3 ) Q %} , (2.17) 

where we have directed the slow variable A* along the 3-d spatial axis such that A^ = CS^i- 
Here the indices a,b = 1,2 are transverse color indices and the index m — 1, 2 is the 
transverse spatial index; X aa = e a/3 A^. The fermion charge of the function (I2.17P is 2. 
For larger N, the fast ground state wave function involves iV(iV — l)/2 such fermion 
factors, each factor carrying the charge 2. The total fermion charge of ipA slov , (^fast) is thus 
-Pfast(-^V) = iV(iV — 1). Speaking of xt^siow), it may carry fermion charges from F s i ow = 
to -F s i ow = 2(N — 1). All together we have 2 2( - N ~^ families of continuum spectrum states 
carrying fermion charges from iV(iV — 1) to (N + 2)(N — 1). Their average charge is 
iV 2 — 1, the half of the maximal fermionic charge -F max - 

The partition function for the system with continuous spectrum is infinite. For exam- 
ple, for the Hamiltonian H = p 2 /2, 

^e-^ 2 = lM , (2.18) 
2ti V 2n 
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with L — > oo. However, the average energy (E)t defined as in (12 .4p does not depend on 
the infinite factor L and is equal to T/2. The Hamiltonian (12.161) involves (D — 1)(N — 1) 
degrees of freedom and we thus obtain 

t n (D-l)(JV-l)/2 

(2.19) 

./V 

(where an infrared regulator \x carrying dimension of energy is introduced) and 

{E)T „ (D-1W-DT p 2Q) 

This result has nothing to do with the supergravity prediction (11.91) ! 

The estimate (I2.19P for the partition function contradicts, however, path integral es- 
timate. The latter is definitely correct (and hence the former is definitely wrong) at high 
temperatures, T ^> A 1 / 3 , when fermions and higher Matsubara modes decouple and the 
partition function is given by the semiclassical estimate (12.121) . the same as for the pure 
YM system [the fermions could only affect the coefficient c in (12.141) ]. 

The paradox is resolved by noting that the spectrum of our system involves besides 
the continuous spectrum also the discrete spectrum with normalized states. Consider first 
4D theory. In this case, the normalized discrete spectrum states of pure bosonic theory 
represent also eigenstates of the full Hamiltonian: the fermion term (11. 7p gives zero when 
acting on the states of zero fermion charge. These normalized eigenstates have the energy 
~ iV 2 A 1//3 as in (12. 9p . The characteristic gap between the lowest and excited normalized 
states in the sector F = is of order A 1 / 3 . Acting on these states by supercharges Q a , we 
can obtain normalized eigenstates of the full Hamiltonian in the sectors F — 1 and F = 2. 
By the same token, one can construct the states in the sectors F = -F max = 2(N 2 — 1), 
F = -F max — 1 and F = F max — 2. Little is known about the structure of normalized 
eigenstates in the sectors with other values of F. It is natural to suggest, however, that 
they also exist and that some of these states (probably, in the sectors with F ~ F max /2 
) may have energy as low as A 1 / 3 (without the N 2 factor). 

If the continuum states did not exist, the partition function at high temperatures 
would be given by the estimate (12.121) (with D — 4) and (E)t by the estimate (I2.13p . 
One can observe now that the latter is much larger than (I2.20p at large N. Heuristically, 
this means that at large N the average energy is determined by the discrete spectrum 
states, while the continuum states are irrelevant. 

Thinking a little bit more in this direction, one could judge that this heuristic impres- 
sion is wrong because one cannot just add the estimates (12.131) and (12.201) and observe 
that (12.131) dominates. One should add the contributions to the partition function rather 
than to the energy. The continuum spectrum contribution to the partition function (12.191) 
involves an infinite factor ^t^- 0-1 )^ -1 )/ 2 and always dominates. 

Thinking still more, one finds, however, that the limits \x — > and iV — > oo do not 
commute. 

• At finite iV and small enough /i, the continuum contribution (12.191) to the partition 
function dominates and one can forget about discrete spectrum. 
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• At finite /i and large enough N, the discrete spectrum contribution (12.121) to the 
partition function dominates and one can forget about continuum. For large enough 
temperatures, T ^> A 1 / 3 , one can also forget about continuum at small N down to 
N = 2. 

Thus, introducing the infrared regulator and playing with this parameter, one can 
get rid of the contribution ( 12.201) . If iV is large enough, one should be able to do it for 
any temperature and not necessary for high temperatures T > A 1//3 where the estimate 
(12.131) is derived. In numerical calculations [7, 8] , no infrared regulator was introduced, 
but the algorithm was chosen such that the continuum spectrum effects were effectively 
suppressed. The functional integral for Z was done by Metropolis algorithm with initial 
values of all components A\ chosen to be of the same order and not very large. It was 
then observed that, for small N, this configuration is unstable such that the field variables 
tend to smear along the flat directions. However, for larger N, the system penetrates the 
valley only after a considerable number of iterations. The larger is N and/or T, the more 
stable is the system. An effective barrier is erected. For large N, one can thus evaluate the 
averages before the system penetrates through this barrier and escapes along the valley. 

Unfortunately, no numerical calculations for the QD system have been done yet, while 
existing calculations for the AD systems [20] are not good enough to conclude about the 
behaviour of (E)t at low T and large N. The authors of [7,8] concentrated on studying 
the 10D system, where they were able to compare their results with the supergravity 
predictions. It seems to us very important to perform the measurements with large enough 
N and good enough statistics also for 4D and 6D systems and compare the results with 
those obtained in the 10.D case. 

What predictions can be made for the behaviour of (E)t in the 4D and 6D cases ? 
Let us assfime[f|the pattern of discrete spectrum states spelled out above: the lowest such 
state has the energy ~ A 1 / 3 and higher excited states behave roughly in the same way as 
in the purely bosonic matrix models up to an overall shift 

E™ ~ N^ 3 — E^ M ~ . (2.21) 

It follows then that, in the limit iV — > oo, with the continuum spectrum effects filtered 
out as explained above, (E)t behaves in the same way as for purely bosonic models, i.e. 
is given by the estimate (12. 13ft at T > A 1 / 3 and approaches zero exponentially fast at 
T — > 0. When iV is large, but finite, this behaviour is valid down to T* ~ A^/lniV 
such that {E)t* ~ NT* and coincides with the continuum estimate (12.201) . At still lower 
temperatures, the system is not contained in the region around A ~ 0, but penetrates 
the barrier and is smeared along the valley. The estimate (12.201) for (E)t is valid in this 
region. 

We are prepared now to go over to 10-D theory and make finally an original remark that 
is raison d'etre for this paper. 10-D theory has one important feature that is lacking in 4D 
and 6D theories: on top of continuum low-energy states and excited discrete spectrum 
states, it involves a normalized vacuum state with zero energy. Its existence was first 

3 We will rediscuss this assumption in Sect. 3. 



8 



discovered in [21] from calculations of integrals for Witten index, but a clearer way to 
see it consists in deforming theory by endowing the scalar fields (in the reduced SQM 
model, they correspond to the components of Aj,...^ with) with mass M [22]. For nonzero 
mass, the classical potential has isolated zero energy minima at large values of the fields. 
Generically, there are several (many) such minima [18], but for the SU(N) groups, there 
is only one classical vacuum. If mass is large, the walls of the potential well around this 
classical vacuum are steep and the wave function of the corresponding quantum vacuum 
is a localized oscillator wave function. By continuity, a localized wave function exists 
also for small values of mass. It is a natural hypothesis that the vacuum state remains 
normalizable also in the limit M — > 0. 

Another argument comes from Born-Oppenheimer analysis of the vacuum wave func- 
tion in the valley. // the normalized state with zero energy exists, its wave function in 
the valley should be represented in the form (12 . 1 5[) with x(x s i ow ) representing the eigen- 
function of the effective Hamiltonian (12.161) with zero eigenvalue. The full wave function 
should be annihilated by supercharges Q a and that means that xi x s\ow) should be annihi- 
lated by effective supercharge acting in Hilbert space of slow variables. One can show that 
normalized solutions to the equation Q^x( x s\ow) = 0, supplemented by the requirement 
of Weyl invariance of x(^siow) following from gauge invariance of the full wave function, do 
not exist in 4D and 6D theories, but the solution exists in the 10D case. It was explicitly 
constructed for SU{2) [23] (see also [18] for pedagogical explanations) and for SU(3) [24]. 
In the simplest N = 2 case, the asymptotic vacuum wave function has the form 

Xvac(A,/i 1 ,..., 8 ) cc (44*™)^™ , (2.22) 

where A = A 3 and (44 ferm )j 7 ' is a fermionic structure representing the 44-plet of SO(9). 

The result (12.221) is obtained in the leading Born-Oppenheimer order. The corrections 
to this result are small when the corrections to the effective Hamiltonian are small. First 
subleading corrections to the effective Hamiltonian are known. Supersymmetry prevents 
the generation of the potential on the valley. In 10D theory, there are also no corrections 
to the metric, i.e. to the term oc E 2 in the Hamiltonian. Q However, there are corrections 
cx E 4 . The exact form of these corrections in the N = 2 case is [27] 

E 2 15 |E| 4 

H efl = 1 - + . . . + terms with fermions. (2.23) 

2 16g 3 |A| 7 

The correction is small iff g|A| 3 ^> 1 (we used E 2 ~ 1/A 2 ). This is also the region 
where the expression (I2.22p for the asymptotic vacuum wave function is valid. On the 
other hand, when g|A| 3 < 1, the separation of slow and fast variables does not work, 
and the wave function depends on all components of A® in a complicated way. Thus, the 
characteristic size of the vacuum wave function is of order A 2 har ~ g~ 2 ^ 3 , which is the 
same as (12.101) . if disregarding iV-dependence there. 

For N > 3, the asymptotic wave function has not been constructed explicitly. We can 
estimate, however, its characteristic size as such A c h ar that the corrections ~ E 4 /A 7 in 

4 This result obtained first in [25] has much in common with AD nonrenormalization theorems [26] . 
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We have 



the effective Hamiltonian are of the same order as the leading term. To this order, the 
effective Hamiltonian is known for any iV [28], 

^w = Ei E "i 2 +iE if A ,_ A L +-. < 2 ^) 

n=l n>m ' 

where we assumed A = diag(A 1 , . . . , A N ) and E = diag(E 1 , . . . , E N ) with J2 n = 
^2 E n = 0. The second term in (I2.24p represents the sum over all positive roots of 
SU(N) (a generalization to an arbitrary group is thus trivial). For large N, it involves 
of order N 2 terms, while the first term has N terms. The estimate for A c h ar is obtained 
from the condition 

N N 2 
~A2 ~ „3 411 • 

^char i/ ^char 

A 2 ^ - ^7F - ^ /9 A- 1/3 • (2-25) 

In Ref. [28], also two-loop corrections to the effective Hamiltonian were evaluated. They 
are estimated as ~ N 3 E S / (g 6 A 1A ) and are of the same order as the leading term at the 
scale (I2.25p . This is probably as well true for higher loop corrections, which should be of 
order 

7?2(n+l) Afn+1 

H n loop ~ . 

We see that the estimated size of the vacuum wave function turns out to be essentially 
larger than the characteristic size (I2.10p of bosonic eigenstates. The large characteristic 
size (I2.25P suggests the existence of the energy scale 

E cZ ~ iV~ 5 / 9 A 1/3 , (2.26) 

which is considerably smaller than the principal energy scale (11.81) . It is natural then to 
assume that the characteristic gap in the spectrum of excited states in 10D SYM quantum 
mechanics is not ~ A 1 / 3 as it was in purely bosonic theory, but smaller, being given by 
the estimate (12.261) . The presence of a large number of discrete spectrum states with the 
energies lying in the interval A 1//3 iV~ 5 / 9 < E < A 1 / 3 modifies essentially the behaviour of 
the partition function. It need not now be exponentially small at T A 1 / 3 , but can well 
display a power behaviour. Assume that 



(Eh ~jv^(j_y 



and determine 7 from the condition {E)t ~ (E con t)T ~ NT at T ~ A^iV -5 / 9 . We obtain 
then 7 = 14/5 in a remarkable agreement with the duality prediction (II. 9p ! 
The law (11.91) implies the behaviour 



T 



9/5 



Z(T) ~ exp { N 2 [ ^ ) -N} (2.27) 
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for the partition function. The normalization (which does not affect {E)t) was chosen such 
that Z(T ~ -E-char) ~ 1- The latter must be true for the discrete spectrum contribution, 
and the continuum spectrum contribution can match it at the scale T ~ i^har ^ choosing 
fi of the same order. When T > E^™, the discrete spectrum contribution dominates. 
Expressing Z{T) into the density of states p(E), 

Z(T) = J p(E)e~ E/T dE , (2.28) 

we see that the behaviour (I2.27j) implies 

p(E)cxexp|iV 2 ^^ ' | . (2.29) 

The critical behaviour (11.91) . (I2.27p . (12.291) should be characteristic in the intermediate 
region 

X 1/3 N- 5/9 <T,E < X 1/3 . (2.30) 
At larger temperatures, the laws (12. 12ft . (12.131) should take over. 



3 Discussion 

Let us summarize our arguments. First, we remark that 10D theory involves besides 
continuum spectrum that is characteristic also for 4D and 6D theories, a normalized 
vacuum state. We estimate a characteristic size of this state by requiring that the loop 
corrections to the effective Hamiltonian at this scale are of the same order as the leading 
term. This gives us a new energy scale (I2.26p . which is lower than the principal energy 
scale (11.81) . Then we conjecture that, on top of the vacuum state, a large family of excited 
states associated with this scale exists. This explains the critical behaviour (11.91) . 

The existence of new scale should also show up in other quantities. For example, 
it suggests that the "extent of space" (12.101) in 10D supersymmetric theory should be 
essentially larger than for purely bosonic system. In particular, the avegage (12.101) should 
grow with N. The existent measurements of ((A®) 2 /N 2 )? in this theory [8] give the 
value that is somewhat larger than in the purely bosonic case, but no growth with iV 
was observed. We do not understand it in the framework of our conjecture and can only 
express a wish that more measurements of this quantity at larger values of N and/or lower 
temperatures were done. Another issue that is not clear now is the range of temperatures 
where the law (11.91) should hold. We have detected only one new scale (I2.26P and no 
other scales. This implies that the law (I1.9P should be valid between these scales, i.e. in 
the range (12.301) . On the other hand, considering this problem on the supergravity side, 
one obtains that subleading in N corrections due to string degrees of freedom become 
essential at T ~ X 1 / 3 ]^- 10 / 21 [7] ) which is somewhat larger than E^ T . It is difficult to 
explain the appearance of this extra scale staying on the matrix model side. 
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Let us go back now to the AD and 6D models. These models do not involve a nor- 
malised vacuum state and we conjectured [see the discussion around Eq. (12.211) ] that the 
pattern of excited normalized states is roughly the same there as in purely bosonic theory. 
On the other hand, there are some indications of the presence of a new energy scale also 
for D = 4,6. When D < 10, the corrections to the moduli space metric do not vanish. 
For an arbitrary gauge group, they have the form [19] 



a D c v 
%|A0')| 3 + 



+ terms with fermions , (3-1) 



where a± = 3/4, a 6 = 1/2 (and a w = 0), J2j * s the sum over all positive roots, = 
a,,(A Cartai1 ), E (j) = aj(E Cartan ) 0, and c v is the adjoint Casimir eigenvalue. For SU(N) 
with large N, the corrections are of order 1 at <?v4 3 har ~ N, which gives 

/V 2 / 3 

D = 4,6: A 2 chai ~ ~ NX- 1 / 3 . (3.2) 

This might be associated with the energy scale ~ A^/iV. //assuming that a family of 
normalised excited states with a characteristic gap ~ X 1 ^ /N is present there, one could 
deduce that the average energy behaves as 

(Ef T DfiD oc T 2 (3.3) 

in the range A 1/3 /iV < T < A 1/3 . 

We would rather lay our own bets not on (13. 3p . but on the scenario spelled out above 
- no low-energy discrete spectrum states and (E)t approaching zero exponentially fast at 
A 1 ^ 3 /(lnA^) < T < A 1 / 3 , with continuum spectrum dependence (E)j> ~ NT taking over 
at still lower temperatures. But theoretical arguments are heuristic and uncertain here 
and only (numerical) experiment can tell us what is true. 

As was mentioned, a nontrivial power behaviour of (E)t is associated with the presence 
of low-energy discrete spectrum states. In principle, one can find these states by solving 
Schrodinger equation numerically. This calculation is, however, much more difficult than 
the calculation of the Euclidean path integral that determines the partition function of 
the system. Up to now, only the simplest AD system with N = 2 was studied [29]. It 
would be very interesting to do it also for higher iV and find out whether a lower energy 
scale ~ \ 1 / 3 /N shows up there. 

A "natural" behaviour of (E)t in quantum mechanics is (E)t oc #d.o.{.T. In our case, 
the power of T is different, which may be associated with the fact, that in the infinite N 
limit, we are dealing actually not with quantum mechanics, but with field theory. Indeed, 
it is known since [30] that the Hamiltonian of the supersymmetric matrix model coincides 
in the infinite N limit with the supermembrane mass operator 



N^L HSQM YM ~ Supermembrane = / d ' V 



(P/) 2 + Xj} 2 + fermionic term 



(3.4) 



3 For example, for SU(N), a„ m (A Cartan ) = A" - A 
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where P[ involves only nonzero modes contribution and {Xi,Xj} = e rs d r Xid s Xj. The 
Hamiltonian (I3.4p is invariant with respect to area-preserving diffeomorphisms (this is 
where gauge symmetry is transformed to in the limit N — »■ oo). Supermembrane theory 
is a (2+l)-dimensional field theory. For the latter, a "natural" law for the area energy 
density is oc T 3 . For sure, this law is derived for a conventional SO (2, 1) invariant field 
theory where, in the limit when interactions are switched off, the spectrum represents a 
tower of oscillators. The model (I3.4|| does not have these features. Still, one can notice 
that 14/5 is numerically close to 3. □. 

We would like to conclude with a general remark. There is a fruitful strategy: when- 
ever you do not understand something in field theory, look at a proper QM system where 
the same phenomenon occurs, analyze it, and you will get chances to improve your under- 
standing. We think that this strategy applies to Maldacena's duality conjecture as well. 
The QM system (11.41) is complicated. Still, it is less complicated than 4D SYM theory at 
strong coupling. If understanding why and how duality works in the former, we will get 
chances to eventually understand it (prove it ) in field theories. 
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Appendix. Supergravity derivation of the law (11.91) . 



To make the paper more self-consistent, we sketch here the derivation of (11.91) on the 
supergravity side. A reader is addressed to the original papers [6] and to the review [32] 
for more details. 

The duality conjecture is that quantum dynamics of different dimensionally reduced 
descendants of 10D SU(N) SYM theory in the large N limit can be accessed by study- 
ing proper classical solutions of 10 D IIA or IIB supergravity. The bosonic part of the 
supergravity action is 



S oc /d iU a;^|e-^[ J R + 4(V0) 2 ]-E^+2| > (A.l) 



where <fi is the dilaton field, F p+ 2 are field strengths of various fundamental (p+ l)-forms 
(Ramond-Ramond forms) that are present in the model, and c p are irrelevant numerical 



6 A further numerological observation is that 5 and 14 are the fourth and the fifth term in the sequence 
of Catalan numbers [31] 
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coefficients. The action (lA.lj) admits black brane solutions, 



ds 



-1/2 



e 
.4 



= f P 

24> K yr&>-3)/2 



-7—p 



n=l 



1/2 



1 - ^— 



7 , + 



•p+i 



irrelevant 



(A.2) 



where / p is a harmonic function in transverse directions (x p +i, . . . , xg). The simplest 
choice is 



fp 



R 7 -p 

r 7-p 



(A.3) 



(r 2 = X^ m =p+i x m anc ^ -R is a constant), (ifig.p is the metric on S 8 ~ p . For example, if 
p = 3 and ro = 0, the metric is reduced to 

,2 



d S 2 



r 



'-dt 2 + dy 2 + dy 2 2 + dyf ) + 



R 2 dr 2 



+ R 2 dQ 



5 ) 



(A.4) 



which is the metric on Ad S§ x S 5 . 

The solutions with p = 3 are relevant when discussing physics of (p + 1) = (3 + 1) - 
dimensional SYM theories. We are interested in the dynamics of (0 + 1) theories and in 
the black hole solution with p — 0, 



r 7/2 



1 " -Z 



+ 



.7/2 



<ir 2 



1-^ 



- + r 2 dnl 



(A.5) 



Let us assume that the black hole size is much less than the characteristic curvature 
radius of the Universe where it sits. This means ro <C R. The black hole HA.5[) has 
a characteristic Hawking temperature and the Bekenstein-Hawking enthropy coinciding 
with the volume of its horizon in Planck units. 

To find the latter, we should simply multiply the factor 



by the factor e 



-2</>(n>) 



oc r, 



21/2 




^(r = r ) oc 



. We obtain 



i 



J/ 2 



S oc r, 



9/2 



(A.6) 

The Hawking temperature is proportional to the so called surface gravity, i.e. gravitational 
acceleration at the horizon (recall the Unruh effect), 

dg 00 



-^Hawking OC ^horizon 



dr 



oc r, 



5/2 



(A.7) 



Combining (I A. 61) and (1A.7|) . we obtain S oc T 9 / 5 and hence (-E)r oc T 14//5 . 



7 Then the metric (|A.2[) is not asymptotically flat and describes a black brane or, which is more 
relevant in the context of etablishing the duality correspondence, a stack of N coinciding black branes in 
the throat. By adding a constant to (|A.3j) . one could obtain an asymptotically flat black brane solution. 
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